Probing the nodal gap in the pressure-induced heavy fermion superconductor CeRhIn 5 
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We report field-orientation specific heat studies of the pressure-induced heavy fermion super- 
conductor CeRhlns. Theses experiments provide the momentum-dependent superconducting gap 
function for the first time in any pressure-induced superconductor. In the coexisting phase of 
superconductivity and antiferromagnetism, field rotation within the Ce-In plane reveals four-fold 
modulation in the density of states, which favors a d-wave order parameter and constrains a theory 
of the interplay between superconductivity and magnetism. 



PACS numbers: 74.20.Rp, 71.27.+a, 74.25. Bt, 75.20.Hr 

The emergence of superconductivity near an antifcr- 
romagnetic quantum critical point (QCP) suggests that 
magnetic fluctuations associated with the QCP may pro- 
vide a glue to form electron pairs In unconven- 
tional superconductors, the superconducting (SC) order 
parameter, which is tied closely to the SC pairing mech- 
anism, breaks the underlying crystalline symmetry and 
contains nodes at which the SC gap becomes zero on 
the Fermi surface. The presence of nodal quasiparticles 
qualitatively changes the density of states (DOS) at the 
Fermi level and leads to a power-law temperature de- 
pendence of thermodynamic properties ■ This power- 
law behavior contrasts with the exponential tempera- 
ture dependence in fully gapped conventional supercon- 
ductors and, therefore, provides an indicator of uncon- 
ventional superconductivity. Unambiguous determina- 
tion of an unconventional SC order parameter, however, 
comes from phase-sensitive and momentum-dependent 
studies that, complemented by the thermodynamic mea- 
surements, have established a d-wave order parameter in 
high- 2^, cuprates Q and p-wave symmetry in Sr2Ru04 
!• 

In heavy fermion superconductors, where the rela- 
tionship between superconductivity and quantum crit- 
icality is most conspicuous [l|, details of the SC gap 
are rarely available because of disorder from chemical 
substitution or pressure environments that are needed 
to induce superconductivity. In this Letter, we re- 
port field-orientation specific heat (C p ) measurements on 
the pressure- induced superconductor CeRhIn 5 . When a 
magnetic field is rotated in the Ce-In plane, four-fold 
modulation in C p is observed deep in the SC state, with 
minima along [100] that indicate gap zeros at this direc- 
tion in a d-wave order parameter. These results provide 
the momentum-dependent SC gap function for the first 
time, which is also the first in any pressure-induced su- 
perconductor, and open a new avenue to insights on the 
pairing mechanism in strongly correlated superconduc- 
tors. 

CeRhlns is an antiferromagnet below Tn = 3.8 K and 
belongs to a family of heavy-fermion compounds CeMIns 



(M=Co, Ir, Rh) that become superconducting for M=Co 
and Ir at 2.3 and 0.4 K, respectively @, @, 0|- Derived 
from the antiferromagnet Celxi3, it crystallizes in the 
tetragonal HoCoGas structure, where a Ce-In plane and 
a M-In block are alternately stacked along the c-axis. 
When subjected to pressure, CeRhIn5 becomes super- 
conducting near 0.5 GPa where Tjv starts to decrease, 
suggesting a correlation between superconductivity and 
magnetism (see Fig. la) [H, @|. As shown in Fig. la, 
there are two critical pressures at 1.75 and 2.3 GPa de- 
noted PI and P2. Below PI, superconductivity and mag- 
netism coexists on a microscopic scale, while a purely SC 
phase only appears above PI 0. P2 is a field-tuned 
quantum critical point hidden by the SC dome Q. At 
pressures above PI, the zero-field nuclear spin-relaxation 
rate 1 jT\ and specific heat show a power-law temperature 
dependence below T c , which is similar to what is found 
in the ambient-pressure heavy-fermion superconductor 
CeCoIns with a d-wave SC order parameter 1(3, 11 1 . In 
the coexisting phase below PI, the temperature depen- 
dence of \jT\ crosses over from T 3 to T- linear with de- 
creasing temperature in the SC state, suggesting a change 
from d-wave singlet pairing to an exotic SC order param- 
eter [T3|. Singular quantum fluctuations in the vicinity 
of an AFM QCP also have been proposed to favor uncon- 
ventional order parameters ranging from a p-wave singlet 
to extended d-wave with additional nodal points 

Whereas power laws in \jT\ and C p are consistent with 
a nodal gap in the coexistence phase, these measurements 
are insensitive to the location of nodes. To resolve this 
important question, we have developed a field-rotation 
specific heat technique under pressure. Single crystals 
used in this study were obtained from a batch where a 
residual resistivity ratio (RRR) was measured routinely 
to be of the order of 500, indicating exceptional sample 
quality. A hybrid clamp-type pressure cell made of Bc- 
Cu/NiCrAl and silicone-fluid transmitting medium were 
used to produce hydrostatic pressure environments up to 
2.5 GPa. The superconducting transition temperature of 
Sn was measured inductively to determine the pressure 
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FIG. 1: (color online) a Phase diagram of CeRhlns as a func- 
tion of pressure [IfJ . Pressure induces a bulk superconducting 
(SC) state at 0.5 GPa that coexists with the antiferromag- 
netic (AFM) phase up to PI, where T c becomes equal to Tn- 
The critical pressure marked by P2 is the field-tuned quan- 
tum critial point hidden by the SC dome, b Heat capacity of 
CeRhlns at 1.47 GPa. The incommensurate AFM transition 
at 2.86 K is followed by a SC transtion at 1.04 K. When 2 T 
is applied (circles), the magnetic transition is barely changed, 
while the superconductivity is completely destroyed. 



at low temperatures. Magnetic field rotation was per- 
formed by using a triple-axis vector magnet with a large 
bore of 2.5" (American Magnetic Inc.) that accommo- 
dated a 3 He cryostat with the pressure cell. Specific heat 
under pressure was measured by an ac calorimetric tech- 
nique where ac heating incurs oscillation in the sample 
temperature (T ac ) that is inversely proportional to the 
sample's heat capacity, i.e., T ac oal/C [l6|. The oscillat- 
ing temperature was measured by a Au-Fe/chromel ther- 
mocouple, which was attached to one face of the plate- 
like sample, while a constantan heater was glued to the 
opposite face. Magnetic field effects on the thermoelec- 
tric voltage of the Au-Fe thermocouple were determined 
independently. 

Figure lb shows the temperature dependence of the 
heat capacity of CeRhlns in its coexisting phase at 
1.47 GPa. Long-range magnetic order with wave vec- 
tor Q=(0.5, 0.5, 0.297) appears at 2.86 K (T N ) and is 
robust against magnetic field; whereas, the SC transition 
occurs at 1.04 K (T c ) and is completely destroyed at 2 T 
(H || a-axis). The magnetic transition at this pressure is 
as sharp as it is at ambient pressure, indicating a negli- 
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FIG. 2: (color online) a Low-temperature heat capacity (C p ) 
of CeRhlns at 1.47 GPa and in low applied fields indicated 
in the legend. Solid and open symbols represent Cp/T for 
field along [100] and [001], respectively, b Temperature de- 
pendence of the upper critical field H C 2- Squares, circles, 
and diamond symbols correspond to fields along [100], [110], 
and [001], respectively. The dotted (or dash-dotted) line is 
a guide for a linear-ff dependence, c Magnetic field depen- 
dence of C/T at 300 mK for H ||[100] at 1.47 GPa (squares) 
and 2.3 GPa(circles) . Here H c2 is 1.2 and 7.7 T for 1.4 and 
2.3 GPa, respectively. The solid line describes least-squares 
fitting by AC/T = C/T - C(0T)/T = AH n with n = 0.68, 
while the dashed line is for n — 0.5. 



gibly small pressure gradient in the pressure cell. 

The low-T heat capacity of CeRhlns is shown in Fig. 2a 
for representative fields along [100] (solid symbols) and 
along [001] (open symbols). These data allow a deter- 
mination of anisotropy of the upper critical field H c2 . 
Within the ab-plane, the upper critical field is es- 
sentially isotropic (see Fig. 2b), but H£ 2 along the c-axis 
is smaller by a factor of two, which reflects the tetragonal 
crystalline structure as found in CeCoIns 



17[ . For mag- 



netic field along any crystalline axis, H C 2 reveals a linear- 
T dependence down to the lowest experimental tempera- 
ture (0.3 K), which is different from CeCoIns where Pauli 
paramagnetic effects dominate conventional orbital fields 
(H£ b = -0.73T C dH c2 /dT) and H c2 is smaller than H% b 
by a factor of 5 [18]. Considering the exceptional sample 
quality of CeRhlns (po = 30 nfl ■ cm), the lack of Pauli 
limiting may arise from the presence of AFM that allows 
a SC order parameter other than a spin singlet [l9j |. We 
also note a positive curvature in near T c , which may 
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be ascribed to multi-band effects since there are three 
major bands that contribute to charge conductivity in 
CeRhIn 5 [Io|. 

Field dependence of the specific heat of CeRhlns does 
not conform to expectations for conventional type-11 su- 
perconductors. Figure 2c shows representative field- 
swept data at 300 mK for H along [100] at 1.47 GPa 
(squares) and at 2.3 GPa (circles). In fully gapped con- 
ventional superconductors, C/T, which is proportional 
to the density of states (DOS), linearly depends on mag- 
netic field because contribution to the DOS comes from 
quasiparticles generated inside vortex cores and the num- 
ber of vortices is proportional to H . In unconventional 
superconductors with nodes, however, C/T increases as 
VH because Doppler-shifted quasiparticles leak through 
the nodes and contribute to the DOS [2l|. A power- 
law form of C/T = C{0T)/T + AH n was used to fit 
the field-dependent specific heat of CeRhIn 5 at 1.47 GPa 
and the solid line in Fig. 2c best describes the data with 
n = 0.68 ± 0.04 for field along [100]. This anomalous be- 
havior may reflect the interplay between superconductiv- 
ity and coexisting antiferromagnetism [TBI, • The field 
dependence of C/T at 2.3 GPa is strikingly different from 
that in the coexisting phase, showing a first-order like 
sharp feature at H C 2 ■ The upper critical field of CeCoIns 
also is first order below T c /2, suggesting that CeCoIris is 



similar to the high-pressure phase of CeRhlns |8j, [23j, [24 1 . 

Figure 3 shows the field-angle specific heat of CeRhIii5 
at 1.47 GPa for field rotations <j> (9) within (perpendicu- 
lar to) the Ce-In plane. In the paramagnetic state at 3 K, 
C/T is featureless as a function of angle (not shown). In 
the AFM state but above T c (squares in Fig. 3a), C/T 
shows a broad feature around [010] that may reflect the 
presence of AFM domains below Tjy (=2.86 K). Upon 
decreasing temperature below T c (=1.04 K), a new mod- 
ulation in C(4>)/T appears in the mixed SC state (circles 
in Fig. 3a). The four- fold modulation is distinctly differ- 
ent from the spectrum observed in the magnetic state, 
indicating that the new structure is not related to mag- 
netism, but arises from the SC gap symmetry. A polar 
sweep at the same pressure, temperature and field (Fig. 
3b) shows only a two-fold modulation that reflects the 
factor of two anisotropy in H C 2 between the ab-plane and 
c-axis. 

Field-rotation specific heat, which is proportional to 
the angle-dependent DOS at low temperature, probes 
the momentum-dependent SC gap function and identi- 
fies gap -zero directions in unconventional superconduc- 
tors [23, 2(1, 27, 2^]. In the Abrikosov state, supercurrent 
circulates around a vortex core and extended quasiparti- 
cles experience a Doppler energy shift, SE « v s • Isf, 
where kp is the Fermi momentum of the quasiparti- 
cles and v s is the supercurrent velocity. Derived from 
the Doppler effects that depend on the relative mag- 
netic field orientation against nodal points on the Fermi 
surface, the DOS precisely reflects the anisotropic gap 
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FIG. 3: (color online) Field-rotation specific heat of CeRhlns 
at 1.47 GPa and 300 mK. a Azimuthal rotation ((f)) of the 
magnetic field (0.5 T) within the Ce-In plane both in the AFM 
(1.8 K) and the coexisting phase of AFM and SC (0.3 K). b 
Polar rotation (8) of the magnetic field (0.5 T) perpendicular 
to the Ce-In plane in the coexisting phase at 0.3 K. The solid 
line is a least-squares fit to a two-fold modulation model: C = 
C +A{l+A 2 \cos(j)\) with C = 265, A = 59.3, and A 2 = 0.30, 
where the oscillation arises from H C 2 anisotropy. 



and modulates with field-rotation angle such that the 
DOS is minimal for field alon g no dal points but is max- 
imal for field along antinodes [25[. Contributions to the 
field-rotation specific heat can be written as Ctotai{4>) — 
Co + Ca ((/>) + Ci((f)), where the constant Cn is the zero- 
field specific heat arising from thermally excited quasi- 
particles and phonons, the two-fold C2{4>) comes from 
sample misalignment, and C±(4>) describes modulation 
from the SC order parameter. The field-angle <f> was de- 
termined against the crystalline a-axis or [100]. The lack 
of a two-fold oscillation indicates that the ab-plane is 
well matched with the field-rotation plane. Figure 4a is 
a plot of the field-induced heat capacity divided by tem- 
perature AC/T at 0.3 K and 1.47 GPa as a function of 
magnetic field angle, where AC = C to tai — Cn = C4 ((/>). 
Minima in AC can be approximated by a parabolic form 
A(l — AiCos4<f)) (dashed line) that is expected for a model 
gap structure of A(9, 0) = An\/1 + cosicfj on a three- 
dimensional Fermi surface. Here 9 is a polar angle deter- 
mined against the c-axis or [001]. 

Figure 4b shows the four-fold oscillation amplitude A\ 
at 0.3 K and 1.47 GPa as a function of magnetic field, 
which is obtained by a least squares fit to the sinusoidal 
form (dashed line in Fig. 4a). Compared to negligible H C 2 
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FIG. 4: (color online) a Field-rotation specific heat of 
CeRhIn 5 at 0.3 K for 0.2 (squares) and 0.9 T (triangles). 
Dashed lines are least-square fits by a sinusoidal form of 
the angular dependence (see text), b Oscillation amplitude 
A4, as a function of the normalized field h(— H/H C 2) with 
H c2 = 1.2 T at 0.3 K. 



anisotropy within the ab-plane (Fig. 2b), the large oscilla- 
tion amplitude (w 3.3%) in AC(cf>) of CeRhlns indicates 
that the four-fold modulation is not due to Fermi surface 
anisotropy, but due to symmetry of the SC gap func- 
tion 



23, Q. Like CeCoIn 5 , minima in AC of CeRhIii5 
are located along < 100 >, indicating d xy symmetry with 
line node along that direction [28|. In a recent theoret- 
ical treatment, however, Vorontsov and Vekhter argued 
that a reversal in the DOS anisotropy occurs at a critical 
field due to additional pair breaking from vortex scat- 
tering; consequently, minima along < 100 > in C p can 
be interpreted instead as evidence for a d x 2_ y 2 gap [2^ |. 
A key future experiment to distinguish the two d-wave 
order parameters is to identify the critical field at which 
anisotropy reversal occurs in the field-rotation specific 
heat of CeRhlns . 

These results constrain possible candidates of the SC 
gap in the coexisting phase of antiferromagnetism and 
superconductivity of CeRhIn 5 . Four- fold modulation of 
the in-plane specific heat excludes an extended d-wave 
with additional nodes at points where the antiferromag- 
netic Brillouin zone crosses the Fermi surface, which was 
suggested to account for the T-linear 1/Ti at lowest tem- 
peratures 14(. A gapless p-wave spin-singlet state that 
was predicted [lj| to be favored over a d-wave symme- 



try near a QCP is also inconsistent with our observations. 
The two-fold modulation in specific heat under polar field 
sweeps, as shown in Fig. 3b, rules out the possibility of 
a hybrid gap with nodes at the poles, which has been 
suggested for Celrln 5 [3l| . Instead, these results indicate 
a d-wave symmetry with line nodes along the c-axis. 

To summarize, four-fold and two-fold modulations in 
representative in-plane and out-of-plane field-rotation 
specific heat measurements on CeRhlns favor a d-wave 
SC order parameter over other suggested candidates in 
the vicinity of the quantum critical point, thus constrain- 
ing a theory of the interplay between superconductivity 
and antiferromagnetism. This work will assure further 
studies of other pressure-induced superconductors where 
the nature of the SC order parameter has yet to be iden- 
tified. 
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